Helium Clustering in Neutron-Rich Be Isotopes 
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Measurements of the helium-cluster breakup and neutron 
removal cross sections for neutron-rich Be isotopes ^^~^^'^'''Be 
are presented. These have been studied in the 30 to 42 MeV/u 
energy range where reaction measurements are proposed to 
be sensitive to the cluster content of the ground-state wave- 
function. These measurements provide a comprehensive sur- 
vey of the decay processes of the Be isotopes by which the 
valence neutrons are removed revealing the underlying a-a 
core-cluster structure. The measurements indicate that clus- 
tering in the Be isotopes remains important up to the drip- line 
nucleus ^''Be and that the dominant hehum-cluster structure 
in the neutron-rich Be isotopes corresponds to a-Xn-a. 

PACS number(s): 25.60. Dz, 25.70.-z, 25.70.Mn, 27.20. -fn 



The structure of nuclear matter at the limits of sta- 
bility has always been a source of fascination. These 
limits are bounded by the extremes of excitation energy 
or spin, on one hand, and by the extremes of either to- 
tal number of nucleons or ratios of protons to neutrons 
(isospin) on the other. Most notably, radioactive beam 
facilities provide access to the limit of neutron excess. 
Here, structural properties are known to differ from nor- 
mal nuclear matter, as they are strongly influenced by 
the weak binding of the valence neutrons. The classic 
example is the halo, where weakly bound neutrons exist 
as a diffuse cloud around a compact core [1] . In this case 
there is a decoupling of the valence neutrons from a core, 
or cluster, of normal nuclear matter. 

Recently, it has been speculated that such clustering 
may be a general feature of drip-line nuclei caused by the 
tendency of the nucleus to distribute the neutron excess 
amoung the core nucleons [2] . In order to maximize the 
overlap of the excess neutrons with protons, the core may 
either deform enhancing the surface area, or optimally, 
form clusters. Thus, the ground states of neutron drip- 



line nuclei may resemble clusters embedded in a cloud of 
weakly bound neutrons. 

Certain nuclei already possess a predisposition towards 
clustering. For example, *Be clearly demonstrates a 
strong a-a cluster character [3,4]. Recent experimen- 
tal and theoretical developments indicate that this un- 
derlying cluster structure has a strong influence on the 
neutron-rich Be isotopes. In particular, the two-centred 
nature of the potential in which the valence neutrons re- 
side may produce molecular type characteristics in these 
nuclei [5,6]. Calculations using the antisymmeterized 
molecular dynamics (AMD) framework suggest such fea- 
tures can occur in the Be isotopes [7,8]. The same model 
predicts that the normally compact boron isotopes be- 
come highly clustered at the drip- line [9] , in line with the 
ideas concerning the optimization of the proton valence- 
neutron overlap [2]. 

Direct access to such structural properties of the 
ground states of very neutron-rich nuclei is extremely 
challenging due to the fact that the ^He cluster breakup 
thresholds increase significantly towards the drip-line. 
However, the AMD model has been used to explore how 
the fragmentation of beams of boron isotopes may be 
used to extract a signature for the clustering in the 
ground state [10]. This study revealed that such cluster- 
ing may be probed via fragmentation if energies in the 
region of 30 MeV/nucleon were used. However, studies of 
the most neutron-rich boron isotope (*^B), as suggested 
in [10], is beyond the capability of present day experi- 
mental facilities. It is, however, possible to study the 
beryllium nuclei all the way from stability to the drip- 
lines. Here we present comprehensive measurements of 
^He-cluster, which were explored theoretically for the 
boron isotopes, and neutron breakup cross sections for 
the beryllium isotopes A=10, 11, 12 and 14. These data 
should provide a significant test of our understanding of 
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the structure of nuclei at the neutron drip-line. 

The measurements were performed at the GANIL ac- 
celerator facility with the production of ^'^'^^^■■^'^■^'^Bc sec- 
ondary beams via the fragmentation of ^'*C and ^'^O pri- 
mary beams. The reaction products were mass, charge 
and momentum analyzed using the LISE3 spectrometer 
producing beams of purity >95% for the lighter isotopes 
and ~20% for the heaviest and count rates of the order 
of 10^ pps for masses 10, 11 and 12 (limited by the count 
rate capacity of the detection system) and ^^50 pps for 
^^Be. Identification of the beam of the beam particle was 
achieved using time-of-flight through the LISE separator 
as measured by a PPAC at the entrance to the reaction 
chamber. The energies of the 4 beams were 30.9, 41.7, 
41.8 and 34.4 MeV/u in order of increasing mass. 

The beam was tracked onto 20 to 275 mg/cm^ car- 
bon targets using two drift chambers, providing a mea- 
surement of the incident position and angle of the pro- 
jectile. The beam and reaction products then entered 
a zero-degree telescope formed from two 500 /um thick, 
16-strip position-sensitive silicon detectors placed 15 cm 
from the target. These two detectors were arranged with 
orthogonal strips, providing a measurement of the in- 
cident ions to <1 mm in both the x and y directions 
(the 0-coordinate being the beam direction). Behind the 
strip detectors was a close packed array of 16 1.5 cm 
thick, 2.5x2.5 cm^, Csl scintillators. These were placed 
at 30 cm from the target so as to cover the same solid 
angle as the strip detectors as measured from the tar- 
get. Calibration of the energy response of the silicon 
and Csl detectors was achieved using a-sources and a 
mixed beam of light-ions of known energy. Neutrons 
produced in reactions of the beam were detected using 
an array of ^100 neutron detectors (DcMoN) arranged 
in a similar manner to that shown in [11]. These de- 
tectors gave a single-neutron detection efficiency of ~15 
% (as in [13]). Monte Carlo simulations of the response 
of the charged particle detection system indicates that 
the efficiency for the detection of the breakup of the Be 
isotopes into two ^Hc nuclei was 40 to 50 % and al- 
most independent of excitation energy of the decaying 
system. The measured relative energy of the decay prod- 
ucts {Erei = M''^rei/2) thus allowcd the excitation energy 
to be deduced {E^ = Erei + Ethresh)- 

A knowledge of the breakup yields, the target thick- 
ness, detection efficiency and the integrated beam ex- 
posure allowed the cross sections for the various breakup 
processes to be calculated. In addition, as noted above, it 
was possible to deduce the excitation energy of the decay- 
ing system. As an example. Fig. 1 shows the ^He-f-^He 
invariant mass spectrum {Ex) for the ^^'^''Be beams. It is 
clear that a large fraction of the events involving the dis- 
sociation of ^'*Be (^^Be) proceed via 6 (4) neutron emis- 
sion to form ^Be, i.e. the decay process does not always 
follow from the decay into neutron-rich helium isotopes 
followed by neutron emission, i.e. unbound excited states 



of ^He. The evidence for this is the observation of peaks 
in the ^Be excitation energy spectrum at energies which 
correspond to the decay of the ground-state (92 keV) and 
the first excited state (3.04 MeV, 2+). There is one fur- 
ther feature in this spectrum, at ~600 keV, which does 
not correspond to a resonance in ^Be but is produced by 
the decay of the 2.43 MeV (5/2") state in ^Bc to the 
tail of the 2+ state in *Be. The presence of this feature 
suggests a fraction of the yield proceeds via sequential 
6n emission to ^Be. Similar features are also present in 
the ''He-l-^He decay spectra measured with the other Be 
projectiles. 

The mass identification provided by the zero degree 
telescope allows the resolution of all the He isotopes and 
thus permits all stages of the decay process to be recon- 
structed. For example, in the case of -'^^Be it is possi- 
ble to reconstruct the decay into ''He+^He, or following 
the emission of neutrons from the projectile, the decay 
of i^Be* into '^Re+'^Re (or ^He+^Hc), or the decay of 
^"Be* to ^He-|-**He, or finally, as discussed above, ®Be 
to *He-|-^He. The cross sections for these processes are 
plotted in Fig. 2. These cross sections are also presented 
together with the ^Be-|-ln (i.e. core fragment plus one 
neutron) cross sections in Table 1. The neutron cross 
sections have been calculated by selecting beam veloc- 
ity neutrons {En >0.37Ebeam/A). In the case of ^^'^^Be 
these were calculated by integrating the neutron angular 
distributions with a Lorentzian single neutron line-shape, 
whilst subtracting the background from reactions within 
the Csl component of the telescopes, as deduced from 
measurements with no target [13]. In the case of ^"'^^Be 
data no background measurements were performed and 
thus the cross sections provides upper limits only. For 
^^Be, measurements of these cross sections, including 
the background subtraction, have already been made by 
Labiche et al. [12,13] at 35 MeV/u, and these are also 
shown in Table 1. Given an analysis of the no target yield 
from the ^^Bc data (~20-30% of events in 2n removal 
channel are from the target) the cross sections are consis- 
tent. Comparison of the cross sections with earlier stud- 
ies show reasonable agreement. The 0.18±0.11 b ^^Be 
In removal cross section compares with that extracted 
by Anne et al. [14] of 0.12±02 b, and the 0.11±0.04 b 
"'^^Be In removal is in close agreement with that measured 
by Navin et al. [15], 0.045±0.05 b, albeit at 78 MeV/u. 
Similar agreement is found with the helium breakup mea- 
surements reported in [16]. 

The measurement of the multi-neutron removal chan- 
nels via the detection of the core plus one neutron 
presents the advantage that the detection efficiency does 
not become prohibitively small. Naively, the true cross 
section for the A — Xn channel should be reduced by 
a factor X. However, the neutron angular distributions 
vary from decay step to decay step and thus the effi- 
ciency is not constant, and, moreover, the multiplicity of 
the emitted neutrons does not necessarily equal the num- 
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bcr of missing neutrons, with the possibihty of projectile 
neutrons interacting strongly with the target [13]. In- 
deed, the measured neutron multiplicities for ^^Be^^^Be 
and i4Be^i"Be are found to be 1.63±0.26 and 2.9±0.8 
[12,13]. Although the error bars are large these multi- 
plicities indicate that there is a tendency for there to be 
slightly less neutrons in the final state than anticipated. 
Nevertheless, for the sake of the comparison with the 
measured helium breakup cross sections the '*~"^Be+ln 
cross sections are plotted in Fig. 2 divided by X, thus 
indicating the strength with which the bound states of 
the A — Xn nucleus are populated. 

It is clear that neutron removal dominates even in 
the case of the nucleus ^"Be in which the neutron decay 
threshold (6.8 MeV) is almost commensurate with that 
for a-decay (7.4 MeV). In all cases the breakup of the 
projectile into helium fragments is several orders of mag- 
nitude smaller than the neutron removal cross sections. 
The helium breakup cross sections steadily increase as 
the neutron removal process feeds states above the clus- 
ter decay thresholds. At the point at which the systems 
arrive at ^Be the cluster breakup and neutron removal 
processes coincide. For each step the ratio of yield in the 
A — Xn and helium breakup gives an idea of the rela- 
tive excitation probabilities to states above the neutron 
and cluster decay thresholds, and subsequent stages re- 
veal the probability that the daughter A — Xn system is 
formed in either bound states or states unbound to either 
neutron of helium decay. These decay systematics should 
allow detailed modeling of the decay process via which 
these nuclei lose all of their neutrons leaving the valence 
a-a cluster core. 

Fig. 3 shows the mass dependence of the cluster decay 
cross sections. It is not possible to plot the data for ^^Be 
in this instance as the main decay channels ^He-|-^He 
^He+^He lead to a *He-|-n-|-^He final state which can- 
not be unambiguously distinguished from ncTitron re- 
moval followed by a-dccay of ^'^Be. Figs 3a and 3b show 
the cross sections for the helium cluster breakup of the 
projectile (first-chance breakup), i.e. ^"'Be^^'He+^He, 
i^Be^^^He-^^He, "'Be^^He+^He, and the total Hc+He 
breakup cross section respectively. The first chance 
breakup cross section is essentially mass independent. 
These systematics indicate that the di-cluster content of 
the nuclei remains relatively constant from ^°Be through 
to ^''Be. However, the total breakup cross section for 
A=10, 11, 12 and 14 shows a strong minimum at A=ll 
increasing to maximum values for '^'^Be and ^''Be. The 
ratio of first chance to total cluster breakup thus reveals 
a maximum at A=12. 

Early calculations of the degree of clustering in the 
beryllium isotopes using the AMD framework indicated 
a decrease in clusterization from ^Be to ^^Be with a 
small increase for ^''Be [8], with something like a fac- 
tor of 2 change predicted in the clusterization between 
the ground states of ^°Be and ^^Be. However, most re- 



cent studies using this framework now indicate low-lying 
states posses strong cluster symmetries [17]. If, as sug- 
gested by the AMD calculations for the reactions with 
the boron isotopes, the breakup process is strongly influ- 
enced by the overlap of the cluster content of the ground 
state and that of excited states above the decay thresh- 
old, then the present measurements indicate that there 
is no change in the degree of clustering in the mass range 
10 to 14. 

By far the most dominant He-cluster breakup channel 
is that corresponding to the removal of all of the valence 
neutrons leaving the a-a core (Fig 2). Table 2 shows 
the energy thresholds for ^'Be+Xn and first-chance He- 
cluster breakup. Remarkably, even though the thresholds 
for removing all valence neutrons are higher in energy, it 
is this channel which dominates, suggesting a strong in- 
fluence of the a-Xn-a structure. Such a result cannot 
be explained by decay phase-space alone or the fact that 
with high numbers of valence neutrons there arc many 
more partitions for the system to decay into. For ex- 
ample, in the case of ^^Be the decay threshold for 6n 
emission lies 4.2 MeV higher than that for cluster decay, 
and inelastic excitation probabilities would be expected 
to decrease with excitation energy. The ^^Be nucleus 
must emit one of six possible neutrons to resonant states 
above the ^^Be 5n decay threshold (10.3 MeV), rather 
than to lower energy states. In turn ^^Be must decay to 
resonances above the ^^Be 4n decay threshold. Thus, al- 
though there is a factor of 6 enhancement in the number 
of possible decay paths for the neutron emission, there 
are 5 intermediate steps in which the system may decay 
to states below the a+a+xn threshold. Moreover, the 
phase space at each step is highly constrained. 

It is possible that the excitation mechanisms leading to 
the first-chance He-cluster breakup and neutron removal 
are radically different. Certainly, direct elastic breakup 
and absorptive processes dominate in the first step over 
resonant breakup as they have been used extensively to 
determine the ground state structure of neutron-rich nu- 
clei [1]. It is thus probable that a significant fraction 
of the a-|-a-|-6n yield for ^^Be is produced via such pro- 
cesses. 

The measurements presented here provide evidence for 
the existence of di-cluster structures in ^""^^^^^Be. More- 
over, they appear to demonstrate that these nuclei pos- 
sess a strong structural overlap with an a-Xn-a configu- 
ration. The comprehensive measurements of the neutron- 
removal and cluster breakup for the first time provide 
experimental data whereby the structure of the most 
neutron-rich Be isotopes can be modeled via their re- 
actions. 
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2.97 1.12 4.29 4.79 11.60 13.27 
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6.81 8.48 


"Hc+^Hc 9.09 
■^Hc+^He 8.95 
"'Hc+SHe 8.81 
*He+6He 7.41 



Figure cind Table Captions 

Table 1 . Cross sections measured for He-cluster breakup 
and core plus one neutron in units of mb. Note in the 
case of ^°'^''Be the neutron cross sections provide upper 
limits. The row marked ^"'Bet are the background sub- 
tracted cross sections from [12,13]. 

Table 2. Energy thresholds for cluster and neutron de- 
cay in units of 
MeV. 

Figure 1. The excitation energy spectrum correspond- 
ing to a+a coincidences from the ^^Be beam. The peaks 
correspond to the the decay of the ground-state (92 keV) 
and the first excited state (3.04 MeV, 2+). The peak 
at ~600 keV is produced by the decay of the 2.43 MeV 
(5/2") state in ^Be to the tail of the 2+ state in ^Be. The 
inset is the excitation energy spectrum corresponding to 
a+a coincidences from the ^"^Be beam. 
Figure 2. Normalized neutron-removal (open dia- 
monds) and cluster breakup (filled circles) cross sections 
for the neutron-rich Be isotopes. For ^^Be the filled dia- 
monds correspond to the measurements of [12,13]. Note, 
the points for the ^Bc breakup (open circles) arc lower 
limits deduced from the 600 keV peak in the ^Be spec- 
trum in Figure 1. The vertical dotted lines indicate the 
mass of the projectile. 

Figure 3. (a) The first-chance cluster breakup cross 
sections, i.e. I'^Be^^He-F^Hc, ^^Bc^ (f'Hc+f'He and 
^Hc+SHe) and ^"Bc^'^Re+^Re, (b) the total He-hHe 
breakup cross sections. 
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